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ABSTRACT 

We present a spectroscopic study of the shape of the dust-extinction law between 1.0 
and 2.2 /mi towards a set of nine ultracompact Hn regions with Ay ;> 15 mag. We 
find some evidence that the reddening curve may tend to flatten at higher extinctions, 
but just over half of the sample has extinction consistent with or close to the average 
for the interstellar medium. There is no evidence of extinction curves significantly 
steeper than the standard law, even where water ice is present. Comparing the results 
to the predictions of a simple extinction model, we suggest that a standard extinction 
law implies a robust upper limit to the grain-size distribution at around 0.1 — 0.3 /im. 
Flatter curves are most likely due to changes in this upper limit, although the effects 
of flattening due to unresolved clumpy extinction cannot be ruled out. 
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1 INTRODUCTION 

Accurate knowledge of the extinction law is crucial to most 
forms of observational astronomy. A detailed understanding 
of extinction in the infrared (IR) waveband, as it applies in 
differing environments, is especially important for the study 
of any heavily embedded object, from young stars to the 
cores of starburst galaxies and AGN. The average extinc- 
tion law in the diffuse interstellar medium is well defined 
in the near-IR and apparently shows little variation (e.g. 
Mathis 1990). However, the broad studies which produce 
this conclusion are based on low-extinction data (e.g. Sav- 
age & Mathis 1979 with A v < 2; Cardelli, Clayton & Mathis 
1989 with Ay < 5). There is a little accurate information 
on how extinction curves might vary along lines of sight to 
heavily embedded objects. It is this latter question that we 
are investigating in this paper. 

There are many reasons to expect an altered extinction 
curve in regions of high gas density and at large Av. Average 
grain sizes are expected to be higher in sheltered molecular 
cloud environments than in the diffuse interstellar medium 
because of the processes of grain coagulation and the growth 
of icy mantles where Av > 3 (Whittet et al. 1983, Tielens 
1989). There is now plenty of evidence for significantly al- 
tered grain-size distributions and increased average grain 
sizes in dense molecular clouds and star-formation regions. 
Such alterations cause variations in the extinction param- 
eter R = Av/E(B — V) and the wavelength of maximum 
polarisation by dichroic absorption where the extinction is 
higher (Savage & Mathis 1979, Whittet et al. 2001). In gen- 



eral, the value of R rises from the interstellar average of 3.1 
to between 4 and 5 in dark clouds and regions of star forma- 
tion. Such a change indicates a flattening of the extinction 
curve at wavelengths shorter than 0.5 /im and is associated 
with alterations in the population of smaller dust grains (e.g. 
Cardelli & Clayton 1991) . 

Additional evidence for grain processing in dense re- 
gions comes from, e.g., IR reflection nebulae associated with 
star formation, in the wavelength dependence of the scat- 
tered intensity (Castelaz et al. 1985), the scattering optical 
depth (Yamashita et al. 1989) and the shape of the 3.09 /im 
ice absorption feature (Pendleton, Tielens & Werner 1990). 
All appear to be affected by the presence of significant num- 
bers of large dust grains, and depletion of smaller grains 
may cause changes in the near-IR continuum extinction in 
low-mass YSOs (Moore & Emerson 1994). 

We have carried out observations of hydrogen recombi- 
nation lines in a number of compact and ultracompact Hn 
regions, and calculated the extinction law for these objects 
by comparing the line ratios with the predictions of well- 
established theoretical models (Hummer & Storey 1995). 
This method has been used previously by Landini et al. 
(1984) and Lumsden & Puxley (1996). The latter found that 
the G45.12+0.13 H 11 region displays an extinction law some- 
what steeper than the standard curve. This paper extends 
the latter study, looking for evidence of widespread vari- 
ations from standard extinction towards ultracompact Hn 
regions with Av J> 15. 
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Table 1. Targets, observed positions and slit orientations 



Object 


IRAS name 


RA (J2000) 


Dec. (J2000) 


Distance 11 (kpc) 


slit angle 6 


G29.96-0.02 


18434-0242 


18 46 03.9 


-02 39 21.9 


6.0 


60° 


G35.20-l.74 


18592+0108 


19 01 46.5 


+01 13 24.8 


3.0 


90° 


G43.89-0.78 


19120+0917 


19 14 26.2 


+09 22 33.9 


9.0 


90° 


G45.45+0.06 


19120+1103 


19 14 21.4 


+11 09 14.1 


4.2/8.8 


90° 


W51d 


19213+1424 


19 23 39.9 


+14 31 08.3 


7.0 


121° 


ROW 57 


11097-6102 


11 11 54.0 


-61 18 23.9 


3.1 


90° 


G298.23-0.33 


12073-6233 


12 10 00.9 


-62 49 54.5 


11.0 


90° 


G307.57-0.62 


13291-6249 


13 32 31.2 


-63 05 19.8 


2.7 


90° 


G337.95-0.48 


16374-4701 


16 41 08.1 


-47 06 46.7 


2.9 


90° 


He2-117 




13 42 36.3 


-61 22 26.6 




90° 


Hc2-96 




15 05 59.5 


-55 59 18.5 




90° 



"See §3.2 for references; 'East of north 



2 OBSERVATIONS 

Observations were made on 18, 19 and 20 May 1997 and 10 
and 12 June 1998 with the common-user infrared imaging 
spectrometer IRIS on the AAT with the cs/36 secondary. 
Nine Hll regions were observed with two planetary nebulae 
as comparison sources. Standards were observed at similar 
airmass, and where possible G-type stars were used, as they 
have few intrinsic absorption features. Where F-type stars 
had to be used, we interpolated over absorption lines in their 
spectra. 

Observations were made in two modes with IRIS. 
The main data presented here were obtained using cross- 
dispersed echelle grisms to obtain complete coverage of IJ 
or HK in a single exposure. The slit width was 1.5" and the 
slit length 13.6", with a pixel scale of 0.8" per pixel. The 
resolution averaged over the wavelength coverage of each 
grism is approximately 400. Sky subtraction was achieved 
through observations of an offset blank sky field after every 
target observation, with the exceptions of the planetary neb- 
ulae, G35.20-l.74 and G43.89-0.78, where the objects were 
compact enough that it was possible to nod along the slit. 
Objects were aligned on the slit using the imaging mode. 
The co-ordinates given in Table 1 represent the mid-point 
of the object as placed in the slit. 

The data were flat-fielded using a tungsten lamp before 
the appropriate on-off pairs were combined. Corrections for 
the curvature of the echelle orders were then applied. Finally 
spectra were extracted in a region covering the peak of the 
emission and approximately 5 rows either side, or in the 
cases where sky subtraction was derived from nodding along 
the slit, spectra were extracted at the peak of the positive 
and negative images by summing over 3 pixel rows and the 
two spectra were then subtracted. Wavelength calibration 
was achieved by means of an argon arc lamp. 

The other observing mode used only the single-order 
H-grism. This gave a resolution of approximately 100 with 
the cs/36 secondary, and wavelength coverage spanning 1.2— 
2.1/nm. These data were acquired to assist in aligning the 
separate IJ and HK echelle spectra, since it was found 
that the poor atmospheric transmission in the region 1.35- 
1.45/nm, caused by water and CO2, made it difficult to join 
the IJ orders to the HK orders accurately. The low reso- 
lution spectra by comparison give a well defined continuum 
shape covering wavelengths included by both echelles. These 
spectra were all acquired during the 1998 observations. The 



Table 2. Hi recombination-line line fluxes for the two planetary 
nebula targets, as a percentage relative to F(Br7) 



Line 


A(air)//im 


Hc2-96 


He2-117 


Case B 


Pa<5 


1.005 


81.5 ± 6.0 


56.7 + 3.9 


201.1 


Pa7 


1.095 


158 ± 11 


121.2 + 8.5 


327.6 


Pa/3 


1.282 


369 ± 26 


329 ± 24 


589.1 


Brl6 


1.556 


5.8 + 0.4 


5.2 + 0.4 


8.0 


Brl5 


1.570 


7.1 ± 0.5 


6.9 + 0.5 


9.7 


Brl4 


1.588 


9.3 + 0.7 


9.0 + 0.6 


11.9 


Brl3 


1.611 


12.1 ± 0.9 


10.9 + 0.7 


14.9 


Brll 


1.681 


19.6 + 1.4 


18.5 + 1.3 


24.7 


BrlO 


1.736 


29.4 + 2.1 


27.1 ± 1.9 


33.1 


Br7 


2.166 


100 + 5 


100 + 5 


100.0 



same nod patterns were used here as in the echelle obser- 
vations. Wavelength calibration was by means of identifying 
the bright Paschen and Brackett lines in the planetary neb- 
ulae. 

The reduction of the straightened IRIS frames then 
followed a standard procedure, using the FIGARO data re- 
duction package. The same procedure was applied to the 
standard-star spectra. Each object spectrum was then di- 
vided by its corresponding standard spectrum, and corrected 
for the black-body shape of the standard. 

The numbers of lines detected depends on the bright- 
ness of the target, the extinction towards it and the atmo- 
spheric seeing. Fluxes for the Hi emission lines were ob- 
tained by fitting Gaussian functions. Errors are taken from 
the fit parameters or assumed to be 5%, reflecting the es- 
timated accuracy of continuum matching between spectral 
bands, whichever is the larger. 



3 RESULTS 

Figure 1 shows the reduced JHK spectrum of G298.23-0.33 
as an example of the data obtained. The separated contin- 
uum level has been matched using the low-resolution grism 
spectrum as described above. 

Hydrogen recombination-line flux ratios are given, rela- 
tive to that of Br7, in Tables 2 and 3. Table 2 contains results 
for the two planetary nebulae, observed as low-extinction 
comparisons to the Hll regions. These were acquired to test 
the accuracy of our observations since they have well de- 
fined optical extinction values. Table 3 contains the main 
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Figure 1. The reduced 1.0 - 2.8-jiim spectrum of G298.23— 0.33 containing lines of the Paschen and Brackett series 
of hydrogen. The three sections have been matched using a low-resolution grism spectrum as a template (see §2). 



Table 3. Measured hydrogen recombination-line fluxes for the Hll regions, as a percentage relative to F(Br7). 



Line 


G29.96-0.02 


G35.20-l.74 


G43.89-0.78 


G45.45±0.06 


W51d 


RCW57 


G298.23-0.33 


G307.57-0.62 


G337.95-0.48 


Pa<5 












14.0 ± 1.0 


3.1 ±0.2 






Pa7 


1.03 ±0.07 








3.77 ±0.26 


40.2 ±2.8 


10.7 ±0.7 


9.5 ±0.7 


2.7 ±0.2 


Pa/3 


16.3 ± 1.1 


14.6 ±1.0 


14.2 ± 1.6 


11.4 ±1.3 


32.9 ± 2.3 


127 ±9 


58.1 ±4.1 


53.9 ±3.8 


20.6 ± 1.4 


Brl6 


1.02 ±0.07 


1.02 ±0.07 






1.92 ± 0.13 










Brl5 


1.30 ±0.10 


1.18 ±0.08 






2.43 ± 0.17 




3.8 ±0.3 






Brl4 


2.10 ±0.15 


1.75 ±0.12 






3.87 ±0.27 


6.36 ± 0.83 


3.9 ±0.5 






Brl3 


2.81 ±0.19 


2.49 ±0.17 






3.85 ± 0.27 


6.8 ± 1.0 


4.6 ±0.3 






Brll 


7.78 ± 0.53 


5.57 ±0.39 


5.45 ± 0.38 


6.46 ± 0.45 


8.53 ± 0.60 


13.7 ± 1.1 


10.9 ±0.8 


10.1 ±0.8 


7.4 ±0.5 


BrlO 


13.06 ±0.88 




6.57 ±0.46 


7.00 ±0.50 


11.6 ± 0.8 


20.0 ± 1.4 


14.5 ±1.0 


17.9 ± 1.3 


9.5 ± 0.7 


Br7 


100 ±5 


100 ±5 


100 ±5 


100 ±5 


100 ±5 


100 ±5 


100 ±5 


100 ±5 


100 ±5 



results for the Hll regions. The results are shown graphi- 
cally in Figures 2 and 3 for the Hll regions and planetary 
nebulae, respectively. The Brl2 line at 1.641 /im has been 
excluded from all fits and calculations since it is blended in 
these observations with the [Fell] 1.644 /im line. 

Observed line ratios are compared to predicted values 
from the models of Hummer & Storey (1995) with n e = 
10 4 cm~ 3 and T e = 10 4 K. Following Landini et al. (1984), if 
the wavelength dependence of the dust opacity is assumed to 
follow a power law, so that k(A) oc A~™, then the differential 
optical depth between two lines, measured as 



Ari; 



In 



[J(A2)//(Ai)] ob 



[J(A 2 )//(Ai)] theory ' 
is equivalent to 



Ar 12 = r Al [(Aa/Ai)-" - l] , 



and the parameters t\ 1 and a can be fitted to the data, in 
this case, by minimising \ 2 using the Marquardt method 
(Fortran routine mrqmin from Press et al. 1993). We have 
taken Ai = A(Br7) = 2.166 /im as the reference wavelength 
in all cases. The results of these fits are listed in Table 4 and 
plotted in Figures 2 and 3. 



3.1 General Results 

Figure 4 shows the compiled results from fitting power-law 
extinction curves to the data. Also plotted is the result for 
G45.12±0.13 obtained by Lumsden & Puxley (1996). Four 
objects have a values more than 2 a below the mean inter- 



4 T. J. T. Moore et al. 



u 

pq 

I 



4 - 



2 - 



- 



4 - 



2 - 



- 



2 - 



- 




1 1.5 2 

wavelength (/xm) 



Figure 2. Extinction-law fits to the HH-region recombination-linc data. 



Table 4. Extinction-law fit parameters 
Source a r(Br7) 



Y 2 



no" 



G29. 96-0.02 


1.90 


± 


0.05 


6.16 


± 


0.07 


4.71 


1.2 


G35.20-l.74 


1.11 


± 


0.12 


4.68 


± 


0.23 


0.48 


4.6 


G43.89-0.78 


1.63 


± 


0.14 


3.32 


± 


0.21 


6.5 


2.8 


G45.45+0.06 


1.33 


± 


0.16 


3.89 


± 


0.26 


10 


2.7 


W51d 


1.95 


± 


0.07 


1.59 


± 


0.07 


2.57 


1.7 


RCW 57 


1.39 


± 


0.14 


1.38 


± 


0.13 


0.57 


1.7 


G298.23-0.33 


1.59 


± 


0.08 


1.76 


± 


0.09 


1.86 


1.7 


G307.57-0.62 


1.79 


± 


0.14 


1.50 


± 


0.13 


1.3 


2.3 


G337.95-0.48 


1.23 


± 


0.09 


3.67 


± 


0.20 


1.9 


2.5 


He2-96 


1.93 


± 


0.32 


0.27 


± 


0.06 


0.49 




He2-117 


1.71 


± 


0.23 


0.46 


± 


0.07 


0.53 





^See §4 



stellar value of 1.8 (e.g. Martin & Whittet 1990) and three 
of these are the three objects with largest fitted values of 
r(Br7). None of the current targets has a significantly larger 
than 1.8. 

The nine Hn regions have 10-/im silicate optical depths 
(tio) estimated by Simpson & Rubin (1990: hereafter SR) 
or Faison et al. (1998), from fits to low-resolution spectra. 
This gives us an independent check of the apparent trend in 
Figure 4. Figure 5(a) shows a plotted against these, taking 
the lower of the no values obtained by SR. Figure 5(b) shows 
our fitted values of r(Br7) plotted against no- 

Five objects have rio values and independent Ay esti- 
mates available in the literature. These suggest that Ay ~ 
IOtio and so no should be comparable in size to the optical 
depth at the wavelength of Br7. 
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Figure 3. Extinction-law fits to data for the PN He2-96 and 
He2-117. The numerical results are given in Table 2. 



3.2 Notes on individual objects 

G29.96-0.02 This object has a likely distance of 6kpc 
(Pratap, Megeath & Bergin 1999). Peeters et al. (2002) re- 
port no detection of 3.1-/iin water-ice absorption in ISO 
data. Estimates of extinction are Ak = 2. 14 ±0.08 (Watson 
et al. 1997), Ak = 1.6 and Ah = 3.6 (Morisset et al. 2002). 
Herter et al. (1981) estimate no = 1.9 ± 0.4. Faison et al. 
(1998) modelled r(2.2) = 1.9 and no = 1.2. Our fit pro- 
duces a r(Br7) consistent with most of these, and a about 
2cr steeper than the standard value of 1.8. 

G35.20~l.74 At a distance of 3kpc (Chan, Henning & 
Schreyer 1996), this source has Ay > 22 mag (Woodward, 
Heifer & Pipher 1985). Our fitted r(Br7) is about twice the 
value extrapolated from this limit using a standard extinc- 
tion law, and is therefore also consistent with the very flat 
curve we derive (5-6 a below the mean interstellar slope). 
Faison et al. (1998) model r(2.2) = 4.1 and no = 4.6. The 
3.1-/im water-ice feature has been detected (paper in prepa- 
ration). 

G43.89-0.78 The estimated distance for this object is 9 kpc 
and no = 2.8 (SR). Our fitted extinction curve is consistent 
with the standard law. 

G45.45+0.06 D = 4.2 or 8.8 kpc (Chan et al. 1996) and 
no = 2.7 (SR). The mean Ak — 2.5 across the Hn region, 
with a peak Ak = 3.6 (with corresponding Av values of 23 
and 33: Feldt et al. 1998). Our -(Brry) is consistent with the 
latter and a is about 3 a natter than the mean extinction 
law. 
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Figure 4. Extinction-law fit parameters a (the power-law expo- 
nent) against the continuum optical depth at Br7. The dashed line 
represents a = 1.8 the mean exponent for the general interstellar 
medium (Martin & Whittet 1990). The open point represents the 
result for G45.12+0.13 found by Lumsdcn & Puxlcy (1996). 



W51d The most luminous component of the giant Hn re- 
gion W51, at a distance of 7.0 kpc (Takahashi et al. 2000). 
Av = 24 ±3 was determined by Goldader & Wynn- Williams 
(1994). Faison et al. (1998) modelled r(2.2) = 2.6 and 
no = 1.7. Our value for a is 2 a larger than the standard 
law. 

RCW57 (G291.3-0.7) from the catalogue of Rodgers, 
Campbell & Whiteoak (1960). Its distance is given as 3.1 kpc 
(SR). Persson, Frogel & Aaronson (1976) determined Ay = 
15 and SR found a silicate optical depth of no = 1-7 (but 
see §4). Our fit produces a value of r(Br7) consistent with 
these and a rather flat a (3 a below 1.8). Pandey & Ogura 
(1998) find R = 4.4 towards the NGC3603 cluster, indicat- 
ing that the extinction curve is also flatter than the average 
at visible wavelengths. 

G298.23-0.33 This object has estimated distance ~ 11 kpc 
(Chan et al. 1996), A v ~ 14 (Persson et al. 1976) and A K = 
0.8 (Martm-Hernandez et al. 2002). no = 1.7 (SR) but H 2 
ice is not detected (Peeters et al. 2002). Our fitted a is flatter 
than the interstellar mean by about 2a. r(Br7) is about 
twice the above Ak value, but is more consistent with Av 
and no- 

G307.57-0.62 The distance estimate is 2.7 kpc (Chan et 
al. 1996) for this source and no = 2.3 (SR). Our fitted 
extinction law is consistent with the standard interstellar 
value. 
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Figure 5. (a) fitted values of a (the extinction-law exponent) 
against the optical depth in the silicate feature no listed by Simp- 
son & Rubin (1990) and Faison ct al. (1998); (b) fitted r(Br7) 
values against no- Uncertainties in no are taken to be 0.1 in each 
case. 




wavelength (/mi) 



Figure 6. A simple two-component model to demonstrate the 
apparent flattening of the measured extinction law caused by un- 
resolved optical depth structure. Here, there are two sources in 
the beam with optical depth ratios of 1.0 (solid line), 3 (short 
dashed) and 10 (long dashed). 



G337.95-0.48 SR give D = 2.9 kpc and no = 2.5 for this 
object. Our spectra were taken during a period of very poor 
seeing which has resulted in only a few of the Brackett se- 
ries lines being measured. The estimate of a we obtain is 
significantly (6 a) below the canonical value. 

Planetary Nebulae He2-96, He2-117 The derived ex- 
tinction law for these objects is consistent with the mean 
ISM law, although the error bars on our fits are large. These 
objects have measured Balmer decrements of ~ 2.0 (He2-96) 
and 2.7 (He2-117) (Tylenda et al. 1992). Assuming the stan- 
dard extinction law, these values correspond to r(Br7) = 
0.33 and 0.44, respectively, consistent with our fitted values 
(Table 4). 



4 DISCUSSION 

Figure 4 suggests an upper boundary to a and an in- 
creased probability that objects with higher absolute ex- 
tinction (r(Br7) > 2.5 or Ay > 25) will have 1/xm - 2/im 
extinction curves that are flatter than the mean interstellar 
law. The data (including the PN and G45.12+0.13 points) 
have a Spearman rank correlation coefficient of r s = —0.636 
which has an associated probability of arising at random of 



2.6 x 10~ 2 for N = 12. This result lends support to the ten- 
tative report of flattened extinction for Ay > 15 by Kenyon, 
Lada & Barsony (1998), from a statistical analysis of photo- 
metric data in p Oph. It also supports the finding of Racca 
et al. (2002) that flatter near-IR extinction laws are found 
where higher extinction regions are present. 

Figure 5(a) shows an apparently similar distribution of 
a with 10— /jm silicate-absorption optical depth obtained by 
fits to IRAS LRS data (SR) and to low-resolution (R = 60) 
spectroscopy (Faison et al. 1998). Figure 5(b) shows that 
our fitted optical depths vary more or less as expected with 
no, although the significance of the correlation is not high 
(the probability of it arising at random is 3.8 per cent). A 
least-squares fit gives a gradient of 1.0 ± 0.3 and intercept 
0.3 ± 0.7. The large scatter should not be too surprising, 
since there may be large uncertainties in no in these low- 
resolution spectra, due to PAH emission features on either 
side of it. There may also be silicate emission components 
filling in the absorption feature (Aitken 1996) and significant 
positional and beam-size discrepancies with our near-IR ob- 
servations. 

Figure 4 also shows that one target (RCW 57) has a 
low fitted value of a without having a particularly high ex- 
tinction. It should be noted, however, that there are large 
amounts of extended, structured emission (e.g. 2MASS K s 
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images) in this target and that it contains a cluster of ob- 
jects, some of which have localised high extinction (e.g. 
no = 3.7 or Av = 59, Barbosa et al. 2003) close by our 
observed position. 

There are two main potential causes of flatter observed 
extinction curves. One is unresolved optical depth structure, 
the other is a physical change in the dust grain properties, 
especially the size distribution, in the material around the 
source. These are discussed below. 



4.1 Structure within the beam 

Figure 6 shows the results of a simple two-component model 
to demonstrate the effects of extinction structure within 
the beam. Three cases are shown in which the wavelength- 
dependence of the extinction is r(A) cx A~ 1,8 for both com- 
ponents, but where optical depths at 2.166 /im are in the 
ratios 1.0, 3.0 and fO.O. ft can be seen that the apparent 
extinction curve becomes flatter with increasing contrast in 
the unresolved optical-depth structure, as well as departing 
from a single power law. Single power-law fits to the lat- 
ter two cases would produce exponents close to f .5 and 1.2, 
respectively. The cause of this artificial flattening is the pen- 
etration of short-wavelength emission through the regions of 
low extinction. 

As mentioned above, at least one Hn region with a low 
fitted a value (RCW 57) is associated with large amounts of 
extended structure. However, there are two reasons to think 
that the effects of unresolved structure may not be the prin- 
cipal cause of the flatter extinction curves. The first comes 
from the model, which shows that the extinction contrast 
has to be large to produce curves as flat as some of those ob- 
served. Secondly, we might expect that more distant sources 
have more unresolved structure and therefore tend to have 
flatter extinction curves, but we find no correlation between 
fitted a values and source distance. We should, therefore 
examine the effect of variable dust-grain properties. 
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Figure 7. Results of simple Mie-scattering model of the near-IR 
extinction opacity due to spherical grains with a size distribution 
n(a) oc A~ 3,5 . (a) shows the effect of increasing the maximum 
grain size; (b) shows the that there is little effect in altering the 
minimum grain size. The dashed line represents a power law of 
A -1,8 . All models curves are normalised to 1 at A = 1.0 fjxa. 



4.2 Grain properties 

A simple model of extinction by spherical grains with a size 
distribution similar to that deduced by Mathis, Rumpl and 
Nordsieck (1977; MRN) can be used to place rough con- 
traints on the grain-size distribution in these regions. We 
have calculated the opacity of an equal-abundance mix of 
silicate and graphite grains using the Mie scattering code 
bhmie by Bohren & Huffman (1983) and the optical con- 
stants used by Weingartner & Draine (2001). We find that 
changes in the population of smaller grains, as manifested 
in a variable minimum grain radius a m i n to the MRN-like 
distribution (n(a) oc a~ 3 ' 5 ), does not affect the slope of the 
near-IR extinction curve between l^im and 2/xm (Figure 7). 
The shape of the curve is much more sensitive to the maxi- 
mum grain radius a max and becomes flatter both when a max 
is small (because the absence of larger grains means that 
absorption becomes a more important component of the ex- 
tinction) and when a ma x is large (when scattering by the 
larger grains is no longer in the Rayleigh regime). Larger 
grains are the more likely case here, since there is other 
evidence of grain growth in star-formation regions (see In- 
troduction). 



These simple models suggest that it is not easy to obtain 
a near-IR extinction law that is significantly steeper than 
the empirical power law (r oc A~™, a ~ 1.8). A steeper 
slope would require a large proportion of weakly absorbing 
grains in the mix. Thick mantles of pure water ice (depth 
> 0.1 pirn) would have this effect, creating a steep extinction 
curve out to the edges of the strong 3.1-^im feature (ignoring 
the weak features at ~1.5/im and ~2/jm). Water ice exists 
toward one source in our list (G35.20-l.74, section 3.2) but, 
given that steep curves are not found, the mantles must be 
relatively thin (significantly less than O.l^im thick). 

The model results also imply that reproduction of the 
standard extinction-law exponent, especially at wavelengths 
between 1 fim and 2 /xm, requires values of a max to be be- 
tween 0.1 /jm and 0.3 /jm (Fig. 7). Physically, this limit may 
arise from the balance between grain coagulation and de- 
struction processes. Departures from it are more likely where 
there is high line-of-sight extinction, but this is not the only 
variable (local density and temperature should be involved) . 

There was little effect on the model predictions from 
introducing increases in the silicate fraction, the addition 
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of amorphous carbon or fluffy composite grains, or from 
changes in the slope of the size distribution. 

Our results are quite consistent with the predictions of 
the model of Kriigel & Siebenmorgen (1994). These authors 
also predict no significant change of slope of the extinction 
curve between wavelengths of 1 and 2pm for grains with a 
size distribution upper limit of 0.3/im, applicable to dust in 
a dense molecular cloud, (see their Figure 12). An observ- 
able flattening of the extinction law is predicted by their 
models for grain distributions with upper limits > 0.3/im. 
Other studies that fit grain models to normal extinction also 
predict upper limits to the grain radius in a narrow range 
close to 0.2 fj,m (e.g. Weingartner & Draine 2001). 



5 CONCLUSIONS 

Determinations of the extinction law between 1/im and 2/im 
have been made towards nine ultracompact Hn regions with 
Av > 15, using Hi line ratios. Four objects have power- 
law fits r oc \~ a with fitted values of a more than 2 a lower 
than the average interstellar value of 1.8. None of the targets 
shows a significantly steeper extinction law. There is an ap- 
parent trend for the sources with highest extinction to have 
flatter extinction curves but low values of a are found at all 
extinctions. Although unresolved clumpy extinction struc- 
ture can cause apparent flattening of the extinction law, the 
observed shallow curves are likely to be due to changes in 
the dust properties. We will examine this trend further in 
Paper II where we study the extinction law from 2-5/itm in 
a sample of HII regions with average extinction higher than 
those presented here. A simple scattering model shows that 
changes to the maximum grain radius in high column den- 
sity regions are the most likely cause and that the canonical 
a = 1.8 exponent suggests that the maximum grain size is 
tightly constrained to be close to 0.2 /im. This upper limit 
must be the result of balance between grain destruction and 
coagulation processes. Any significant change to the maxi- 
mum grain radius due to a change in this balance will cause 
a flattening of the near-IR extinction curve. 
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